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Abstract 

A model of a shear band as a zero-thickness nonlinear interface is proposed and tested using 
finite element simulations. An imperfection approach is used in this model where a shear band, 
that is assumed to lie in a ductile matrix material (obeying von Mises plasticity with linear 
hardening), is present from the beginning of loading and is considered to be a zone in which 
yielding occurs before the rest of the matrix. This approach is contrasted with a perturbative 
approach, developed for a J 2 -deformation theory material, in which the shear band is modelled 
to emerge at a certain stage of a uniform deformation. Both approaches concur in showing that 
the shear bands (differently from cracks) propagate rectilinearly under shear loading and that 
a strong stress concentration should be expected to be present at the tip of the shear band, two 
key features in the understanding of failure mechanisms of ductile materials. 

Keywords: 


1 Introduction 

When a ductile material is brought to an extreme strain state through a uniform loading process, 
the deformation may start to localize into thin and planar bands, often arranged in regular lattice 
patterns. This phenomenon is quite common and occurs in many materials over a broad range of 
scales: from the kilometric scale in the earth crust (Kirby, 1985), down to the nanoscale in metallic 
glass (Yang, 2005), see the examples reported in Fig. [lj 

After localization, unloading typically^] occurs in the material outside the bands, while strain 
quickly evolves inside, possibly leading to final fracture (as in the examples shown in Fig. [2j where 
the crack lattice is the signature of the initial shear band network) or to a progressive accumulation 
of deformation bands (as for instance in the case of the drinking straws, or of the iron meteorite, or 
of the uPVC sample shown in Fig.[lJ or in the well-known case of granular materials, where fracture 
is usually absent and localization bands are made up of material at a different relative density, Gajo 
et al. 2004). 

It follows from the above discussion that as strain localization represents a prelude to failure of 
ductile materials, its mechanical understanding paves the way to the innovative use of materials in 
extreme mechanical conditions. Although shear bands have been the subject of an intense research 

1 For granular materials, there are cases in which unloading occurs inside the shear band, as shown by Gajo et al. 
( 2004 ). 
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Figure 1: Examples of strain localization. From left to right, starting from the upper part: A merlon in the 

Finale Emilia castle failed (during the Emilia earthquake on May 20, 2012) in compression with a typical ‘X-shaped’ 
deformation band pattern (bricks are to be understood here as the microstructure of a composite material). A 
sedimentary rock with the signature of an ‘X-shaped’ localization band (infiltrated with a different mineral after 
formation). A stone axe from a British Island (Museum of Edinburgh) evidencing two parallel localization bands and 
another at a different orientation. A runestone (Museum of Edinburgh) with several localized deformation bands, 
forming angles of approximatively 45° between each other. A polished and etched section of an iron meteorite 
showing several alternate bands of kamacite and taenite. Deformation bands in a strip of unplasticized poly(vinyl 
chloride) (uPVC) pulled in tension and eventually evolving into a necking. An initially regular hexagonal disposition 
of drinking straws subject to uniform uniaxial strain has evolved into an ‘X-shaped’ localization pattern. A fracture 
prevails on a regularly distributed network of cracks in a vault of the Amiens dome. ‘X-shaped’ localization bands in 
a kaolin sample subject to vertical compression and lateral confining pressure. A thin, isolated localization band in 
a sedimentary layered rock (Silurian formation near Aberystwyth). 
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Figure 2: Regular patterns of localized cracks as the signature of strain localization lattices. From left to right: 
Dried mud; Lava cracked during solidification (near Amboy crater); Bark of a maritime pine (Pinus pinaster); Cracks 
in a detail of a painting by J. Provost (‘Saint Jean-Baptiste’, Valenciennes, Musee des Beaux Arts). 


effort over the last twenty years (see the review given by Bigoni, 2012), many fundamental questions 
still remain open: i.) Why are shear bands a preferred mode of failure for ductile materials? ii.) 
Why do shear bands propagate rectilinearly under mode II, while cracks do not? iii.) how does a 
shear band interact with a crack or with a rigid inclusion? iv.) Does a stress concentration exist at 
a shear band tip? v.) How does a shear band behave under dynamic conditions? 

The only systematic^ attempt to solve these problems seems to have been a series of works by 
Bigoni and co-workers, based on the perturbative approach to shear bands (Argani et al. 2014; 
2013; Bigoni and Capuani, 2002; 2005; Piccolroaz et al. 2006). In fact problems (i.), (ii.), and (iv.) 
were addressed in (Bigoni and Dal Corso, 2008 and Dal Corso and Bigoni, 2010), problem (iii.) in 
(Dal Corso et al. 2008; Bigoni et al. 2008; Dal Corso and Bigoni, 2009), and (v.) in (Bigoni and 
Capuani, 2005). 

The purpose of the present article is to present a model of a shear band as a zero-thickness 
interface and to rigorously motivate this as the asymptotic behaviour of a thin layer of material, 
which is extremely compliant in shear (Section [2]). Once the shear band model has been developed, 
it is used (in Section [3]) to (i) demonstrate that a shear band grows rectilinearly under mode II 
remote loading in a material deformed near to failure and (ii.) estimate the stress concentration 
at the shear band tip. In particular, a pre-existing shear band is considered to he in a matrix 
as a thin zone of material with properties identical to the matrix, but lower yield stress. This is 
an imperfection, which remains neutral until the yield is reached in the shear band. The present 
model is based on an imperfection approach and shares similarities to that pursued by Abeyaratne 
and Triantafyllidis (1981) and Hutchinson and Tvergaard (1981), so that it is essentially different 
from a perturbative approach, in which the perturbation is imposed at a certain stage of a uniform 
deformation process. To highlight the differences and the analogies between the two approaches, 

2 Special problems of shear band propagation in geological materials have been addressed by Puzrin and Ger¬ 
manovich (2005) and Rice (1973). 
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the incremental strain field induced by the emergence of a shear band of finite length (modelled as a 
sliding surface) is determined for a ^-deformation theory material and compared with finite element 
simulations in which the shear band is modelled as a zero-thickness layer of compliant material. 


2 Asymptotic model for a thin layer of highly compliant material 
embedded in a solid 


A shear band, inside a solid block of dimension , is modeled as a thin layer of material (of semi¬ 
thickness /i, with h/H <C 1) yielding at a uniaxial stress ay\ which is lower than that of the 
surrounding matrix ay n \ Fig. [3 Except for the yield stress, the material inside and outside the 
layer is described by the same elastoplastic model, namely, a von Mises plasticity with associated 
flow rule and linear hardening, defined through the elastic constants, denoted by the Young modulus 
E and Poisson’s ratio z/, and the plastic modulus E p: see Fig. Ilf- 



Figure 3: (a) A shear band inside a ductile material modeled as a thin layer of highly compliant 
material ( E ep /E <C 1) embedded in a material block characterized by a dimension i7, such that 
h/H <C 1; both materials obey the same von Mises plasticity model represented by the uniaxial 
stress behaviour reported in (b), but having a different yield stress (lower inside than outside the 
shear band). 

At the initial yielding, the material inside the layer (characterized by a low hardening modulus 
E ep = EE p /(E + E p )) is much more compliant than the material outside (characterized by an elastic 
isotropic response E). 

For h/H <C 1, the transmission conditions across the layer imply the continuity of the tractions, 
t = [t2i,£22] T , w hich can be expressed in the asymptotic form 

|[*2i] = 0(h), [* 22 ! = 0(h), (1) 
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where [•] denotes the jump operator. The jump in displacements, [ix] = [[xq], |^2]] T , across the 
layer is related to the tractions at its boundaries through the asymptotic relations (Mishuris et ah, 
2013; Sonato et ah, 2015) 


*2l([^ll, I'M) = 


V^HP + 4 HP + 6hajf_ -gH| 

(3 E + 2(1 + u)E p )y/Z\u i] 2 + 4H 2 2 h 


+ 0(h), 


( 2 ) 


*22 (HP HI) = 


(E + 2(1 - v)E p y 3HF+4HP + 8h(l - 204 S) E iM 


+ 0(h ), (3) 


(1 - 2 i /)(3 E + 2(1 + u)E p )y/3lu i] 2 + 4[u 2 P 2h 

involving the semi-thickness h of the shear band, which enters the formulation as a constitutive 
parameter for the zero-thickness interface model and introduces a length scale. Note that, by ne¬ 
glecting the remainder 0(h), Eqs. © and © define nonlinear relationships between tractions and 
jump in displacements. 

The time derivative of Eqs. © and © yields the following asymptotic relation between incre¬ 
mental quantities 

1k_ 1 + k 0 (M,M) H, (4) 

where the two stiffness matrices K-\ and Kq are given by 


K i = 


K n = 


E 


2(3 E 2(1 v)E v 
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E-n 


Y 


0 

E + 2(1 — v)E p 
1 - 2v 


HI 2 


-HIM 

hi ] 2 


(5) 


( 6 ) 


(3 E + 2(1 + u)E p )( 3H1 2 + 4H] 2 )3/2 L-H1M 
Assuming now a perfectly plastic behaviour, E p = 0, in the limit h/H —> 0 the condition 

HI = 0 (7) 

is obtained, so that the incremental transmission conditions © can be approximated to the leading 
order as 

t-ljSHH. (8) 

Therefore, when the material inside the layer is close to the perfect plasticity condition, the incre¬ 
mental conditions assume the limit value 

*21 = 0, [* 2 ] = 0, (9) 

which, together with the incremental version of eq. ©2, namely, 

M = 0, (10) 

correspond to the incremental boundary conditions proposed in Bigoni and Dal Corso (2008) to 
define a pre-existing shear band of null thickness. 

The limit relations ^ and (10) motivate the use of the imperfect interface approach (Mishuris, 
2004; Mishuris and Ochsner, 2005; 2007; Antipov et al. 2001; Bigoni et ah, 1998) for the modelling 
of shear band growth in a ductile material. A computational model, in which the shear bands are 
modelled as interfaces, is presented in the next section. 
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3 Numerical simulations 


Two-dimensional plane-strain finite element simulations are presented to show the effectiveness of 
the above-described asymptotic model for a thin and highly compliant layer in modelling a shear 
band embedded in a ductile material. Specifically, we will show that the model predicts rectilinear 
propagation of a shear band under simple shear boundary conditions and it allows the investigation 
of the stress concentration at the shear band tip. 

The geometry and material properties of the model are shown in Fig.[4j where a rectangular block 
of edges H and L > H is subject to boundary conditions consistent with a simple shear deformation, 
so that the lower edge of the square domain is clamped, the vertical displacements are constrained 
along the vertical edges and along the upper edge, where a constant horizontal displacement u\ is 
prescribed. The domain is made of a ductile material and contains a thin ( h/H <C 1) and highly 
compliant ( E ep /E <C 1) layer of length H/2 and thickness 2 h = 0.01 mm, which models a shear 
band. The material constitutive behaviour is described by an elastoplastic model based on linear 
isotropic elasticity (E = 200000 MPa, ^ = 0.3) and von Mises plasticity with linear hardening (the 
plastic modulus is denoted by E p ). The uniaxial yield stress for the matrix material is equal 
to 500 MPa, whereas the layer is characterized by a lower yield stress, namely, — 400 MPa. 


ui 



Figure 4: Geometry of the model, material properties and boundary conditions (that would corre¬ 
spond to a simple shear deformation in the absence of the shear band). The horizontal displacement 
u\ is prescribed at the upper edge of the domain. 

The layer remains neutral until yielding, but, starting from that stress level, it becomes a 
material inhomogeneity, being more compliant (because its response is characterized by E ep ) than 
the matrix (still in the elastic regime and thus characterized by E). The layer can be representative 
of a pre-existing shear band and can be treated with the zero-thickness interface model, Eqs. <© 
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and ([ 3 ]). This zero-thickness interface was implemented in the ABAQUS finite element softwar^] 
through cohesive elements, equipped with the traction-separation laws, Eqs. 0 and Q, by means 
of the user subroutine UMAT. An interface, embedded into the cohesive elements, is characterized 
by two dimensions: a geometrical and a constitutive thickness. The latter, 2/i, exactly corresponds 
to the constitutive thickness involved in the model for the interface 0 and ([5]), while the former, 
denoted by 2 h g , is related to the mesh dimension in a way that the results become independent 
of this parameter, in the sense that a mesh refinement yields results converging to a well-defined 
solution. 

We consider two situations. In the first, we assume that the plastic modulus is E p — 150000 
MPa (both inside and outside the shear band), so that the material is in a state far from a shear 
band instability (represented by loss of ellipticity of the tangent constitutive operator, occurring 
at E p = 0) when at yield. In the second, we assume that the material is prone to a shear band 
instability, though still in the elliptic regime, so that E p (both inside and outside the shear band) is 
selected to be ‘sufficiently small’, namely, E p — 300 MPa. The pre-existing shear band is therefore 
employed as an imperfection triggering shear strain localization when the material is still inside the 
region, but close to the boundary, of ellipticity. 

3.1 Description of the numerical model 

With reference to a square block (L = H = 10 mm) containing a pre-existing shear band with 
constitutive thickness h — 0.005 mm, three different meshes were used, differing in the geometrical 
thickness of the interface representing the pre-existing shear band (see Fig. [ 5 ] where the shear band 
is highlighted with a black line), namely, h g = {0.05; 0.005; 0.0005} mm corresponding to coarse, 
fine, and ultra-fine meshes. 

(a) 


L. 


Figure 5: The three meshes used in the analysis to simulate a shear band (highlighted in black) 
in a square solid block (L = H = 10 mm). The shear band is represented in the three cases as 
an interface with the same constitutive thickness h — 0.005 mm, but with decreasing geometric 
thickness h g \ (a) coarse mesh (1918 nodes, 1874 elements, h g = 0.05 mm); (b) fine mesh (32079 
nodes, 31973 elements, h g = 0.005 mm); (c) ultra-fine mesh (1488156 nodes, 1487866 elements, h g 
= 0.0005 mm) 

The three meshes were generated automatically using the mesh generator available in ABAQUS. 
In order to have increasing mesh refinement from the exterior (upper and lower parts) to the 

3 ABAQUS Standard Ver. 6.13 has been used, available on the AMD Opteron cluster Stimulus at UniTN. 
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interior (central part) of the domain, where the shear band is located, and to ensure the appropriate 
element shape and size according to the geometrical thickness 2 h gi the domain was partitioned into 
rectangular subdomains with increasing mesh seeding from the exterior to the interior. Afterwards, 
the meshes were generated by employing a free meshing technique with quadrilateral elements and 
the advancing front algorithm. 

The interface that models the shear band is discretized using 4-node two-dimensional cohesive el¬ 
ements (COH2D4), while the matrix material is modelled using 4-node bilinear, reduced integration 
with hourglass control (CPE4R). 

It is important to note that the constitutive thickness used for traction-separation response is 
always equal to the actual size of the shear band h — 0.005 mm, whereas the geometric thickness h g: 
defining the height of the cohesive elements, is different for the three different meshes. Consequently, 
all the three meshes used in the simulations correspond to the same problem in terms of both material 
properties and geometrical dimensions (although the geometric size of the interface is different), so 
that the results have to be, and indeed will be shown to be, mesh independent. 


3.2 Numerical results 


Results (obtained using the fine mesh, Fig.^) in terms of the shear stress component a ±2 at different 
stages of a deformation process for the boundary value problem sketched in Fig. [4] are reported in 
Figs. [6] and [7| 
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Figure 6: Contour plots of the shear stress a 12 for the case of material far from shear band instability 
( E p = 150000 MPa). The grey region corresponds to the material at yielding a \2 > 500/a/3 ^2 288.68 
MPa. Three different stages of deformation are shown, corresponding to a prescribed displacement at 
the upper edge of the square domain u\ = 0.037418 mm (a), u\ — 0.037518 mm (b), u\ = 0.037522 
mm (c). The displacements in the figures are amplified by a deformation scale factor of 25 and the 
percentages refer to the final displacement. 

In particular, Fig. prefers to a matrix with high plastic modulus, E p = 150000 MPa, so that the 
material is far from the shear band formation threshold. The upper limit of the contour levels was 
set to the value g \2 — 500/a/3 ^2 288.68 MPa, corresponding to the yielding stress of the matrix 
material. As a result, the grey zone in the figure represents the material at yielding, whereas the 
material outside the grey zone is still in the elastic regime. Three stages of deformation are shown, 
corresponding to: the initial yielding of the matrix material (left), the yielding zone occupying 
approximately one half of the space between the shear band tip and the right edge of the domain 
(centre), and the yielding completely linking the tip of the shear band to the boundary (right). Note 
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Figure 7: Contour plots of the shear stress a 12 for the case of material close to shear band instability 
(. E p — 300 MPa). The grey region corresponds to the material at yielding au > 500/Three 
different stages of deformation are shown, corresponding to a prescribed displacement at the upper 
edge of the square domain u\ — 0.0340 mm (a), u\ — 0.0351 mm (b), u\ — 0.03623 mm (c). The 
displacements in the figures are amplified by a deformation scale factor of 27. 


that the shear band, playing the role of a material imperfection, produces a stress concentration at 
its tip. However, the region of high stress level rapidly grows and diffuses in the whole domain. At 
the final stage, shown in Fig. almost all the matrix material is close to yielding. 

Fig. 0 refers to a matrix with low plastic modulus, E p = 300 MPa, so that the material is close 
(but still in the elliptic regime) to the shear band formation threshold (E p = 0). Three stages of 
deformation are shown, from the condition of initial yielding of the matrix material near the shear 
band tip (left), to an intermediate condition (centre), and finally to the complete yielding of a narrow 
zone connecting the shear band tip to the right boundary (right). In this case, where the material 
is prone to shear band localization, the zone of high stress level departs from the shear band tip 
and propagates towards the right. This propagation occurs in a highly concentrated narrow layer, 
rectilinear, and parallel to the pre-existing shear band. At the final stage of deformation, shown in 
Fig. 03, the layer of localized shear has reached the boundary of the block. 

Results in terms of the shear strain component 712 , for both cases of material far from, and 
close to shear band instability are reported in Figs. [ 8 ] and [9j respectively. In particular, Fig. [ 8 ] 
shows contour plots of the shear deformation 712 for the case of a material far from the shear band 
instability {E p = 150000 MPa) at the same three stages of deformation as those reported in Fig. [ 6 j 
Although the tip of the shear band acts as a strain raiser, the contour plots show that the level of 
shear deformation is high and remains diffused in the whole domain. 

Fig- H shows contour plots of the shear deformation 712 for the case of a material close to the 
shear band instability (E p = 300 MPa), at the same three stages of deformation as those reported in 
Fig- 0 It is noted that the shear deformation is localized along a rectilinear path ahead of the shear 
band tip, confirming results that will be reported later with the perturbation approach (Section [4]). 

The shear deformation 712 and the shear stress a 12 along the x-axis containing the pre-existing 
shear band for the case of a material close to strain localization, E p — 300 MPa, are shown in 
Fig. 10, upper and lower parts, respectively. Results are reported for the three meshes, coarse, fine 
and ultra-fine (Fig. [ 5 ]) and at the same three stages of deformation as those shown in Figs. [T] and 
[9| The results appear to be mesh independent, meaning that the solution converges as the mesh is 
more and more refined. 

The deformation process reported in Figs. 00 and [TO] can be described as follows. After an 
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Figure 8 : Contour plots of the shear deformation 712 for the case of material far from shear band 
instability (E p = 150000 MPa). Three different stages of deformation are shown, corresponding 
to a prescribed displacement at the upper edge of the square domain u\ — 0.037418 mm (a), 
u\ = 0.037518 mm (b), u\ = 0.037522 mm (c). The displacements in the figures are amplified by a 
deformation scale factor of 25. 



Figure 9: Contour plots of the shear deformation 712 for the case of material close to shear band 
instability (E p = 300 MPa). Three different stages of deformation are shown, corresponding to a 
prescribed displacement at the upper edge of the square domain u\ = 0.0340 mm (a), u\ — 0.0351 
mm (b), u\ = 0.03623 mm (c). The displacements in the figures are amplified by a deformation 
scale factor of 27. 


initial homogeneous elastic deformation (not shown in the figure), in which the shear band remains 
neutral (since it shares the same elastic properties with the matrix material), the stress level reaches 
(J 12 = 400/\/3 — 230.9 MPa, corresponding to the yielding of the material inside the shear band. 
Starting from this point, the pre-existing shear band is activated, which is confirmed by a high shear 
deformation 712 and a stress level above the yield stress inside the layer, —5 mm < x < 0 (left part 
of Fig. 10). The activated shear band induces a strain localization and a stress concentration at its 
tip, thus generating a zone of material at yield, which propagates to the right (central part of Fig. 


10) until collapse (right part of Fig. 10) 


In order to appreciate the strain and stress concentration at the shear band tip, a magnification 
of the results shown in Fig. [l0]in the region —0.2 mm < x < 0.2 mm is presented in Fig. mi Due 
to the strong localization produced by the shear band, only the ultra-fine mesh is able to capture 
accurately the strain and stress raising (blue solid curve), whereas the coarse and fine meshes smooth 
over the strain and stress levels (red dotted and green dashed curves, respectively). The necessity of 
a ultra-fine mesh to capture details of the stress/strain fields is well-known in computational fracture 
mechanics, where special elements (quarter-point or extended elements) have been introduced to 
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Figure 10: Shear deformation 712 (upper part) and shear stress 072 (lower part) along the x-axis 
containing the pre-existing shear band for the case of a material close to a shear band instability 
E p — 300 MPa. The black dotted line, in the bottom part of the figure, indicates the yield stress 
level, lower inside the pre-existing shear band than that in the outer domain. Three different stages 
of deformation are shown, corresponding to a prescribed displacement at the upper edge of the 
square domain u\ — 0.0340 mm (left), u\ = 0.0351 mm (center), u\ = 0.03623 mm (right). 


avoid the use of these ultra-fine meshes at corner points. 

For the purpose of a comparison with an independent and fully numerical representation of the 
shear band, a finite element simulation was also been performed, using standard continuum elements 
(CPE4R) instead of cohesive elements (COH2D4) inside the layer. This simulation is important 
to assess the validity of the asymptotic model of the layer presented in Sec. [2| In this simulation, 
reported in Fig. [I2j the layer representing the shear band is a ‘true’ layer of a given and finite 
thickness, thus influencing the results (while these are independent of the geometrical thickness 2 h g 
of the cohesive elements, when the constitutive thickness 2 h is the same). Therefore, only the fine 
mesh, shown in Fig. [ 5 ) 3 , was used, as it corresponds to the correct size of the shear band. The 
coarse mesh (Fig. |5^i) and the ultra-fine mesh (Fig. [5^) would obviously produce different results, 
corresponding respectively to a thicker or thinner layer. Results pertaining to the asymptotic model, 
implemented into the traction-separation law for the cohesive elements COH2D, are also reported in 
the figure (red solid curve) and are spot-on with the results obtained with a fully numerical solution 
employing standard continuum elements CPE4R (blue dashed curve). 

A mesh of the same size as that previously called ‘fine’ was used to perform a simulation of a 
rectangular block (H = 10 mm, L = 4 H = 40 mm) made up of a material close to shear band 
instability ( E p = 300 MPa) and containing a shear band (of length H /2 = 5 mm and constitutive 
thickness 2 h = 0.01 mm). Results are presented in Fig. [l3j In parts (a) and (b) (the latter is a 
detail of part a) of this figure the overall response curve is shown of the block in terms of average 
shear stress 072 — T/L(T denotes the total shear reaction force at the upper edge of the block) and 
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Figure 11 : Shear and stress concentration at the shear band tip. Shear deformation 712 (upper 
part) and shear stress cr 12 (lower part) along the x-axis containing the pre-existing shear band for 
the case of a material close to a shear band instability E p — 300 MPa. Three different stages of 
deformation are shown, corresponding to a prescribed displacement at the upper edge of the square 
domain u\ — 0.0340 mm (left), u\ — 0.0351 mm (center), u\ = 0.03623 mm (right). 


average shear strain 712 = u\/H. In part (c) of the figure contour plots of the shear deformation 
712 are reported at different stages of deformation. It is clear that the deformation is highly focused 
along a rectilinear path emanating from the shear band tip, thus demonstrating the tendency of the 
shear band towards rectilinear propagation under shear loading. 

Finally, the incremental shear strain (divided by the mean incremental shear strain) has been 
reported along the x-axis in Fig. M at the two stages of deformation considered in Fig. [TO] and 
referred there as (a) and (c). These results, which have been obtained with the fine mesh, show that 
a strong incremental strain concentration develops at the shear band tip and becomes qualitatively 
similar to the square-root singularity found in the perturbative approach. 

4 The perturbative vs the imperfection approach 

With the perturbative approach, a perturbing agent acts at a certain stage of uniform strain of 
an infinite body, while the material is subject to a uniform prestress. Here the perturbing agent 
is a pre-existing shear band, modelled as a planar slip surface, emerging at a certain stage of a 
deformation path (Bigoni and Dal Corso, 2008), in contrast with the imperfection approach in 
which the imperfection is present from the beginning of the loading. 

With reference to a x\-X 2 coordinate system (inclined at 45° with respect to the principal pre¬ 
stress axes £/-£//), where the incremental stress tij and incremental strain Sij are defined (i, j =1,2), 
the incremental orthotropic response under plane strain conditions (^3 = 0) for incompressible ma- 
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Figure 12: Results of simulations performed with different idealizations for the shear band: zero¬ 
thickness model (discretized with cohesive elements, COH2D) versus a true layer description (dis¬ 
cretized with CPE4R elements). Shear deformation 712 (upper part) and shear stress 072 (lower 
part) along the horizontal line y — 0 containing the pre-existing shear band for the case of a material 
close to a shear band instability E p = 300 MPa. Three different stages of deformation are shown, 
corresponding to a prescribed displacement at the upper edge of the square domain u\ = 0.0340 
mm (left), u\ — 0.0351 mm (center), u\ — 0.03623 mm (right). 


terials (in +i 22 = 0 ) can be expressed through the following constitutive equations 

tn = 2 /i^n + p, t 2 2 = - 2 /i^n + p, *12 = M*7i2, 


(Bigoni, 2012 J] 
( 11 ) 


where p is the incremental in-plane mean stress, while fi and /i* describe the incremental shear 
stiffness, respectively, parallel and inclined at 45° with respect to prestress axes. 

The assumption of a specific constitutive model leads to the definition of the incremental stiffness 
moduli y and / 1 *. With reference to the ^-deformation theory of plasticity (Bigoni and Dal Corso, 
2008), particularly suited to model the plastic branch of the constitutive response of ductile metals, 
the in-plane deviatoric stress can be written as 

*7 - *// = kei\ei\^ N ~ l \ (12) 


In equation (12) fc represents a stiffness coefficient and N E (0,1] is the strain hardening exponent, 
describing perfect plasticity (null hardening) in the limit N 0 and linear elasticity in the limit 
N —» 1. For the ^-deformation theory, the relation between the two incremental shear stiffness 
moduli can be obtained as 

p* = N/jl, (13) 


4 Note that the notation used here differs from that adopted in (Bigoni and Dal Corso, 2008), where the principal 
axes are denoted by x\ and X 2 and the system inclined at 45° is denoted by x\ and £ 2 . 
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Figure 13: Results for a rectangular domain (L = 40 mm, = 10 mm) of material close to shear 
band instability (E p = 300 MPa) and containing a preexisting shear band (of length H/2 = 5 
mm and constitutive thickness 2 h = 0.01 mm), (a) Overall response curve of the block in terms 
of average shear stress 072 = T/L, where T is the total shear reaction force at the upper edge of 
the block, and average shear strain 712 = u\jH. (b) Magnification of the overall response curve 
012 — 712 around the stress level corresponding to the yielding of the shear band, (c) Contour plots 
of the shear deformation 712 at different stages of deformation, corresponding to the points along 
the overall response curve shown in part (b) of the figure. The deformation is highly focused along a 
rectilinear path emanating from the shear band tip. The displacements in the figures are amplified 
by a deformation scale factor of 50. 


so that a very compliant response under shear (/i* <C /i) is described in the limit of perfect plasticity 

TV 0. 
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Figure 14: The incremental shear strain 712 (divided by the mean incremental shear strain 712 ) 
along the x-axis at the two stages of deformation reported in in Fig. 10 and labeled there as (a) 


and (c). It is clear that a strong strain concentration develops at the tip of the shear band, which 
becomes similar to the square-root singularity that is found with the perturbative approach (Section 
[4] and Fig. [I 6 j 


The perturbative approach (Bigoni and Dal Corso, 2008) can now be exploited to investigate the 
growth of a shear band within a solid. To this purpose, an incremental boundary value problem is 
formulated for an infinite solid, containing a zero-thickness pre-existing shear band of finite length 
21 parallel to the x\ axis (see Fig. 15) and loaded at infinity through a uniform shear deformation 

.coo 

7l2- 



Figure 15: A perturbative approach to shear band growth: a pre-existing shear band, modelled as 
a planar slip surface, acts at a certain stage of uniform deformation of an infinite body obeying the 
J 2 -deformation theory of plasticity 

The incremental boundary conditions introduced by the presence of a pre-existing shear band 
can be described by the following equations: 

t 2 i(a:i,0 ± ) = 0, 1*22(* 1 ,0)] = 0, [^ 2 (^ 1 ,0)] = 0, V|cci| < l. (14) 

A stream function ^(aq,^) is now introduced, automatically satisfying the incompressibility con¬ 
dition and defining the incremental displacements iij as ii\ — ^ 2 , and U 2 = — ^q. The incremental 
boundary value problem is therefore solved as the sum of ^°{x i,aq), solution of the incremental 
homogeneous problem, and ^ p (aq,aq), solution of the incremental perturbed problem. 
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The incremental solution is reported in Fig. [16] for a low hardening exponent, N = 0.01, as a 
contour plot (left) and as a graph (along the aq-axis, right) of the incremental shear deformation 
712 (divided by the applied remote shear qg). Note that, similarly to the crack tip fields in fracture 
mechanics, the incremental stress and deformation display square root singularities at the tips of 
the pre-existing shear band. Evaluation of the solution obtained from the perturbative approach 
analytically confirms the conclusions drawn from the imperfection approach (see the numerical 
simulations reported in Fig. [9] and [l3|) , in particular: 


• It can be noted from Fig. [l6] (left) that the incremental deformation is highly focussed along 
the x\ direction, confirming that the shear band grows rectilinearly; 


• The blow-up of the incremental deformation observed in the numerical simulations near the 
shear band tip (Fig. [14 \ is substantiated by the theoretical square-root singularity found in 
the incremental solution (Fig. [l6j right). 



Figure 16: Incremental shear strain near a shear band obtained through the perturbative approach: 
level sets (left) and behaviour along the aq-axis (right). 


5 Conclusions 

Two models of shear band have been described, one in which the shear band is an imperfection 
embedded in a material and another in which the shear band is a perturbation which emerges 
during a homogeneous deformation of an infinite material. These two models explain how shear 
bands tend towards a rectilinear propagation under continuous shear loading, a feature not observed 
for fracture trajectories in brittle materials. This behaviour is a basic micromechanism of failure 
for ductile materials. These models show also a strong stress concentration at the shear band tip, 
which can strongly promote shear band growth. 
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